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We perform finite-temperature dynamical simulations of the arrest of a rotating Bose-Einstein condensate 
by a fixed trap anisotropy, using a Hamiltonian classical-field method. We consider a quasi-two-dimensional 
condensate containing a single vortex in equilibrium with a rotating thermal cloud. Introducing an elliptical 
deformation of the trapping potential leads to the loss of angular momentum from the system. We identify the 
condensate and the complementary thermal component of the nonequilibrium field, and compare the evolution 
of their angular momenta and angular velocities. By varying the trap anisotropy we alter the relative efficiencies 
of the vortex-cloud and cloud-trap coupling. For strong trap anisotropies the angular momentum of the thermal 
cloud may be entirely depleted before the vortex begins to decay. For weak trap anisotropies, the thermal cloud 
exhibits a long-lived steady state in which it rotates at an intermediate angular velocity. 
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I. INTRODUCTION 

Since the first experimental realizations |Q] |2l [3] of dilute- 
gas Bose-Einstein condensates, there has been much interest 
in the properties of quantum vortices in such systems, and 
the effect of thermal atoms on their dynamics and stability. 
The description of the dynamics of vortices in the presence 
of thermal excitations provides a challenging test for dynami- 
cal theories of cold bosonic gases, and promises new insights 
into the role of thermal excitations in the dynamics of vor- 
tices in systems less amenable to such ab initio descriptions 
00 |6]E]]. The stability of a vortex state in a dilute conden- 
sate was first investigated theoretically by Rokhsar JSJ, who 
showed that a vortex is subject to decay in the presence of 
a nonrotating thermal cloud, providing a physical interpreta- 
tion of a negative-energy excitation of the vortex previously 
found by Dodd et al. 0. Fedichev and Shlyapnikov ifTUl 
then put forward an analytic theory of the dynamics of vortex 
decay in the presence of a nonrotating thermal cloud, based 
on a two-fluid model of superfluidity. Madison et al. ifTTl 
subsequently observed experimentally a nonexponential de- 
cay in vortex survival probability, and the first evidence of an 
increase in the displacement of the vortex from the trap center 
during the decay, which was also observed in the experiments 
of Anderson et al. Ifl2l . Zhuravlev and co-workers lfT3ll pro- 
vided an analytic description of the decay of vortex arrays, 
building on the work of IfTUl . and including the rotational dy- 
namics of a nonstationary thermal cloud on the basis of the 
theory of Guery-Odehn fl4l . Their work predicted two lim- 
iting regimes of relaxation dynamics, depending on the rela- 
tive efficiencies of vortex-cloud and cloud-trap coupling: the 
so-called rotating trap limit, in which the vortex array and 
thermal cloud relax collectively as a single rigid body, and the 
static trap limit, in which the thermal cloud is quickly arrested 
by the trap anisotropy, and nonexponential decay [ 15 1 of the 
array rotation follows. 

Abo-Shaeer et al. fl6| performed experiments in the rotat- 
ing trap regime, and observed the expected exponential decay 
and strong temperature dependence of the decay rate. Rosen- 
busch et al. [ 17 1 observed the decay of a single vortex exper- 



imentally, and found a much less severe dependence on tem- 
perature than that of [16|. They conjectured that the thermal 
cloud was rapidly arrested by residual trap anisotropy, leading 
to the static-cloud vortex-decay scenario of [lOj. Classical- 
field simulations performed by Schmidt et al. 11811 focused on 
the dynamics of a strongly nonequilibrium 'phase-imprinted' 
vortex state. Duine et al. [ 19| presented an analytical descrip- 
tion of the decay of a vortex in a nonrotating thermal cloud, 
obtained from the stochastic field theory of [20 1 using a vari- 
ational ansatz GTl . which was subsequently extended to in- 
clude the effects of cloud rotation by Bradley and Gardiner 

ma. 

In the present paper we consider a condensate containing a 
vortex initially at thermal and rotational equilibrium with a ro- 
tating thermal cloud in a highly oblate trap, which is isotropic 
in the plane. Such a configuration is obtained as an ergodic 
classical-field equilibrium with fixed angular momentum on 
the order of n per atom about the trap axis Il23ll . Due to the 
conservation of angular momentum, this rotating equilibrium 
configuration is stable, provided that the trapping potential re- 
mains invariant under rotations about its axis. We then intro- 
duce an elliptical deformation of the trap which breaks this 
rotational symmetry, leading to the loss of angular momen- 
tum from the atomic field and thus the slowing of the rotating 
cloud, and consequently the decay of the vortex. Our simula- 
tions are the first to describe the arrest of the rotation of both 
the condensed and noncondensed components of the field fol- 
lowing the introduction of a trap anisotropy, and our method 
describes the coupled relaxation dynamics of the two compo- 
nents implicitly. As predicted by [13|, we find that the re- 
sponse of the condensed and noncondensed components of 
the field may be different depending on the relative efficien- 
cies of vortex-cloud and cloud-trap coupling, where the latter 
depends on the ellipticity of the trap deformation. By varying 
the anisotropy over a range of values and monitoring the evo- 
lution of the condensed and noncondensed components of the 
field, we observe a rich phenomenology, ranging from an adi- 
abatic steady state, to violently nonequilibrium dynamics, in 
which the rotation of the thermal cloud essentially decouples 
from that of the condensate. 
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This paper is organized as follows. In Sec. [II] we discuss 
our classical-field formalism, and the parameters of the phys- 
ical system we simulate. In Sec. Ill we discuss the results 



of a representative simulation. In Sec. IV we discuss the de 



pendence of the dynamics of the vortex and thermal field on 
the trap anisotropy. In Sec. [V] we summarize our findings and 
present our conclusions. 



n. FORMALISM 

A. Classical-field method 

The general formalism of (projected) classical-field meth- 
ods has recently been reviewed at length in lt24l . and so 
here we merely outline the particular method we use in this 
work. The system we represent is described by the second- 
quantized, cold-collision (s-wave) Hamiltonian 



U 



j dx*¥ 



t (x)4' t (x)4'(x)4'(x), 

where ^(x) is the bosonic field operator, which satisfies 
[T'(x), 4*' (x')] = 6(x-x'), and U - 4nh 2 a/m is the interaction 
potential written in terms of scattering length a and atomic 
mass m. The single-particle Hamiltonian is of form 



-/z 2 V 2 
2m 



+ V (x) + V e (x) = Hi + V e (x), (2) 



where Vo( x ) = ( OT /2){<w 2 (x 2 + y 2 ) + cj 2 z 2 } is a cylindrically 
symmetric trapping potential and V 6 (x) = £ma> 2 {y 2 - x 2 ) is 
an additional anisotropic potential. In deriving the classical- 
field method, we replace the field operator *T(x) with the pro- 
jected classical field ifr(x) = JineL a n<Pn( x )- The sum here 
is over eigenmodes of the cylindrically symmetric single- 
particle Hamiltonian [H^ p <p n (x) - e„(p„(x)] with eigenvalues 
satisfying e„ < Er, where Er is the cutoff energy. The pro- 
jected classical field is thus a vector in a low-energy subspace 
(L) of the appropriate single-particle Hilbert space. Making 
this replacement we obtain the classical Hamiltonian 



U , 



H CF = \ dx if,*(x)H sp if,(x) + -l-Mx)! 4 . (3) 



■/ 

Defining the projector 

P/W^W f dyf n (y)f(y), 

neL J 



(4) 



we can express the Hamilton's equation for i^(x) obtained 
from Eq. Q as 



Br 



(5) 



which is termed the projected Gross-Pitaevskii equation l24l . 
In the present work we use the rotationally invariant (Gauss- 
Laguerre) projector introduced in ll25l . For a discussion of the 



evaluation of projections of the 2-body interaction and trap 
anisotropy operations in this basis, we refer the reader to ll25ll 
and 11261 respectively. We simply note here that this method 
conserves both the energy and normalization of the field, and 
that the field angular momentum obeys the Ehrenfest relation 
of the continuous field theory 



(6) 



[where the bar denotes spatial averaging: A = 
j dxiff*(x)Aiff(x)], which is a consequence of the relation 
[V, L z ] = 0. We note finally that as our formalism describes 
only the low-energy region L, it does not account for any 
population above the cutoff. The results of any classical- 
field theory are in general cutoff dependent, and the cutoff 
dependence of the formalism used here was quantified in 
ll26l . As noted in that work, ultimately this dependence can 
only be eliminated by the inclusion of above-cutoff physics, 
and currently no prescription for including such effects in 
general nonequilibrium scenarios is known. In practice, 
however, the low-energy region defined by Er ~ 3/u, with fi 
the chemical potential (see Ref. E6ll and Refs. therein) is ex- 
pected to include the most important dynamical mechanisms, 
the above-cutoff atoms influencing the dynamics only at a 
quantitative level. 



B. Simulation procedure 

Following [26 1, we choose physical parameters correspond- 
ing to 23 Na atoms confined in a strongly oblate trap, with 
trapping frequencies (a> r ,a> z ) — 2n x (10,2000) rad/s. The 
radial harmonic confinement defines the units of length (ro = 
■\/h/m(D r ) and time [one trap cycle (eye) = 2nu>~ 1 ] in which 
we will quote our results. The s-wave scattering length is 
a = 2.75n;«, placing our system in the quasi-two-dimensional 
regime (oscillator length L — yjti/ma) z » a), with an ef- 
fective two-dimensional (2D) interaction parameter U2D = 
2 yfhrha / ml z . We are thus justified in representing our system 
by a 2D model, provided that we maintain thermodynamic pa- 
rameters (chemical potential yU and temperature T) such that 
fi + ksT <s ha) z ll27l at all times during the evolution. The 
2D representation is formalized in our classical-field method 
by choosing a cutoff E R = 30fiw,. <K Ticl> z , such that the low- 
energy space L excludes all modes with excitation along the 
z-axis. In all our simulations this energy cutoff (projection 
operation) is effected in an inertial (laboratory) frame, and the 
low-energy space consists of 465 single-particle modes. 

We form a finite-temperature initial state following the pro- 
cedure of ll23l : we construct a (nonequilibrium) randomized 
classical-field configuration over the modes comprising the 
space L, with chosen normalization, energy and angular mo- 
mentum first integrals, and we evolve this state for some 
time (10 4 trap cycles), so that the field has time to migrate 
to an equilibrium configuration. Following [23 1, we take as 
the starting point for forming these configurations the ground 
Gross-Pitaevskii (GP) eigenstate with No = 1.072 x 10 4 
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atoms in a frame rotating at angular velocity Qo = 0.35a> r . 
This state has chemical potential /j g = 10.35/zw,- and energy 
£g = 7.646 x l0 4 fiLo r in an inertial frame, and angular mo- 
mentum Lo = %Nq. To this state we add energy and angular 
momentum, forming a configuration with E = l.lOZig, and 
L = 1.20Lo. We find that the corresponding equilibrium state 
is one in which a single vortex precesses very close to the trap 
axis, i.e., with precession radius r v < 770, where the healing 
length 770 = 0.20ro is estimated from the density of the ground 
state ll23l . In this configuration the vortex is at equilibrium 
with the thermal component of the field, and the angular ve- 
locity Q c of the condensate (which is the same as the cloud's at 
equilibrium) is close to that at the I = (L z )/ttN — > 1~ limit of 
the mechanically unstable ll28l l29l precessing-vortex branch 
QUI . If additional angular momentum is added to this state, 
the condensate mode becomes unstable to the nucleation of a 
second vortex, i.e., the single-vortex state we consider is es- 
sentially saturated with angular momentum for the given val- 
ues of the other conserved integrals (normalization and en- 
ergy). We find for this equilibrium initial state Q c as 0.23w r , 
/i x \2hcj r , k-gj ss 14/zav and condensate fraction f c ss 0.91. 
Having formed this initial state, we evolve field trajectories in 
the laboratory frame in the presence of trap anisotropies in the 
range 0.005 < e < 0.1. This procedure corresponds in each 
case to the sudden introduction of the anisotropy at time t — 0. 
We evolve each of these trajectories with an adaptive integra- 
tor OTI . with accuracy chosen such that the relative change 
in field normalization is < 10~ 9 per time step taken, until the 
vortex is expelled from the condensate and the total angular 
momentum of the field has been lost through its interaction 
with the anisotropy. 



III. RESULTS 



In this section we consider the results of a classical-field 
simulation with a particular choice of the trap anisotropy, 
e = 0.025. The relaxation dynamics of the classical field 
in this case exhibit many features of interest, and serve as a 
point of comparison for simulations with weaker or stronger 
anisotropies. Position-space densities of the classical field at 
representative times are presented in Fig. |TJa-f)- Initially the 
vortex precesses very close to the trap axis [Fig.flla)], and it 
remains close to the trap axis for some time. Fig7]Tjb) shows 
the density at time t = 500 eye, in which the vortex core re- 
mains within ~ 770 of the trap axis. By t m 1000 eye the vor- 
tex has begun to spiral out of the central density bulk which 
contains the condensate. Fig. [Tfc) shows that by t — 1500 
eye the vortex has undergone significant radial displacement. 
This increase in radial displacement continues [Fig.[TJd)], un- 
til the vortex core approaches the violently evolving conden- 
sate boundary [Fig. [TJe)]. At t w 3070 eye the vortex is lost 
into the peripheral thermal material, leaving the condensate 
vortex-free and essentially irrotational [Fig.[TJf)]. 



A. Vortex trajectory 



We now quantify the behavior of the vortex during its de- 
cay. We observe that the vortex executes an essentially spiral- 
like motion as it decays, in agreement with the predictions of 
ifTUll and the simulations of lfl8l l32l . However, the vortex tra- 
jectory here is strongly stochastic, and the vortex orbits the 
trap axis many (~ 700) times during its decay, so we do not 
present the full vortex trajectory. To characterize the radial 
drift of the vortex, we track its location at a frequency of 25 
samples per trap cycle, and average the resulting series over 
intervals of 100 samples (4 eye) in order to smooth out rapid 
fluctuations due to thermal density fluctuations in the back- 
ground field against which the vortex moves and uncertain- 
ties introduced by sampling the vortex location on a Cartesian 
grid. In Fig. |2ja) we present the smoothed vortex radius (r\) 
data. It is clear from Fig. |2ja) that there is an initial lag of 
~ 1000 eye between the introduction of the trap anisotropy at 
t = eye and the beginning of the upward trend in vortex dis- 
placement. We note that the displacement radius does not in- 
crease monotonically, but exhibits large oscillations during its 
increase. The period of these oscillations is often > 100 trap 
cycles, spanning many periods of the vortex orbit, in contrast 
to the trajectories presented in |[T8l [32l . in which the vortex 
generally precesses at a greater radius on each passage about 
the trap center. In Fig.[2|b) we plot the vortex precession fre- 
quency. In practice, to evaluate this frequency, we follow the 
procedure of E3l . forming the covariance matrix (classical 
one-body density matrix) 

p(x,x') = <-AW(x')>n, (7) 

where (■ • ■ )n denotes a time average, in a frame rotating at an- 
gular velocity Q about the trap axis. We vary £2 such that the 
largest eigenvalue of p(x, x') is maximized. The value of Q at 
which the maximum occurs is thus that of the rotating frame 
in which the coherent fraction of the classical field is most 
stationary, which forms a best estimate for the condensate an- 
gular frequency (vortex precession frequency) Q c 11231 . As 
we consider a nonequilibrium scenario, we construct Eq. (j7]) 
by averaging over short time periods [33], in each case calcu- 
lating the average of 250 consecutive classical-field samples 
over a ten cycle period. Fig. [2jb) shows that the vortex pre- 
cession frequency determined in this manner increases as the 
vortex displacement increases, as is well known for a zero- 
temperature condensate [29, 34] and observed in the simula- 
tions of lfT8l l32l . Moreover, we observe that the oscillations 
in vortex radius r v at late times t > 2000 eye are accompanied 
by oscillations in the condensate angular velocity. The oscil- 
lations in the two are positively correlated, i.e., the decreases 
in vortex radius are associated with periods of slowing of the 
condensate rotation during its otherwise steady increase with 
time. 
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B. Rotational dynamics of condensate and thermal cloud 

The above definition of the condensate mode in terms of 
short-time covariance matrix eigenvectors allows us to resolve 
the dynamics of the condensed and noncondensed compo- 
nents of the field. As in 11231 . we introduce the decomposition 
of the classical-field one-body density matrix in terms of its 
eigenvectors \xi) and corresponding eigenvalues n, (indexed 
in order of decreasing eigenvalue), 



P = «oUfO>CkTol + n k\Xk)(Xk\ =Po+ Pth, 



(8) 



/.>o 



which separates it into condensed and noncondensed parts. As 
in ll23l . we then define the averages of a single-body operator 
J in the condensed and noncondensed components of the field 
by </)o = Trfpo-/} and {J) lh = Tr{p th 7) respectively. We use 
this decomposition to calculate the angular momentum of the 
condensate ((L z )o) and thermal cloud ((L ; ) t h). These quan- 
tities are presented in Fig. [3ja), along with the total angu- 
lar momentum of the field ((L z ) tot ). We observe that in the 
first few hundred trap cycles after the introduction of the trap 
anisotropy, the angular momentum of the thermal cloud [solid 
line in Fig. |3ja)] undergoes approximately exponential decay, 
reaching values close to zero by t « 1000 trap cycles. By con- 
trast, the angular momentum of the condensate (dashed line) 
is essentially unchanged during this period, which we iden- 
tify with the initial period of quiescence of the vortex noted 
in Sec. Ill A At t « 1100 eye, the angular momentum of the 
condensate begins to decay, and after t * 3200 eye the con- 



densate is essentially irrotational, with some small value of 
angular momentum remaining in the form of surface excita- 
tions. We note that the angular momentum of the thermal field 
component increases over the period t e [1000, 3200] eye. We 
understand this as follows: Initially the angular momentum of 
the thermal cloud is lost due to its interaction with the trap 
anisotropy. This makes the vortex thermodynamic ally unsta- 
ble, and thus it begins to decay, liberating its angular momen- 
tum to the cloud. The cloud loses its angular momentum to 
the trap at a rate proportional to its angular velocity (frictional 
loss d(L z )/dt oc -Q t h lfl3l ). which is initially small. There- 
fore, the cloud's angular momentum increases as it gains an- 
gular momentum from the vortex faster than it loses it to the 
trap. As the cloud's rotation rate increases the rate of dissipa- 
tion of (L z )± increases, however, the precessing vortex con- 
densate exhibits a thermodynamic anomaly: as the conden- 
sate loses angular momentum, its angular velocity increases 
11281 [29l 30 1 . The rate at which the angular momentum of the 
vortex is dissipated to the thermal cloud is proportional to their 
relative velocities IfTOll , and so in the present case this transfer 
dominates the loss rate of angular momentum of the thermal 
cloud throughout the vortex decay. As shown in Fig. [3jb), 
this leads to 'run-away' spin-up of the vortex precession fre- 
quency, with the thermal-cloud rotation rate being driven up 
in sympathy. 

Closer inspection of the angular momenta (L z )q and (L z ) t h 
reveals negatively correlated oscillations in the angular mo- 
menta of the two components during the decay, as a result 
of the detailed nonlinear dynamics of angular momentum ex- 
change and loss during the decay process. We identify this 
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FIG. 2: (Color online) Quantities characterizing the vortex motion, 
(a) Vortex radial displacement and (b) angular velocity of the con- 
densate orbital. 
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FIG. 3: (Color online) Quantities characterizing the temporal evo- 
lution of the condensate and thermal cloud, (a) Angular momenta of 
the condensate and thermal cloud, and total angular momentum of 
the field, (b) Angular velocities of the condensate and thermal cloud. 



oscillation in angular momentum transfer with the oscillations 
in vortex rotation rate discussed in Sec. Ill A In Fig. [3jb) we 
plot the angular velocity of the condensate Q c (dots) derived 
in the estimation of the condensate fraction (Sec. |III A) , and 
the angular velocity of the thermal cloud Q t h = (-L z )th/(®c)th 
(line), where we assume the expectation value of the classical 
moment of inertia C = r 2 as an estimate of the cloud's true 
moment of inertia. In (23 1 we discussed the level of uncer- 
tainty inherent in this procedure, nevertheless it yields a clear 
qualitative description of the decay dynamics. The oscilla- 
tions in rotation rate of the thermal cloud are clearly visible 
here, and we note that they are positively correlated with the 
oscillations in angular velocity of the condensate, due to its 
anomalous rotational response. We note that oscillatory be- 
havior arises already in a linear analysis of the arrest of a rotat- 
ing Boltzmann gas [ 14 1 by a trap anisotropy. It is therefore not 
surprising that similar oscillations occur in the transfer of an- 
gular momentum from the condensate to the nonequilibrium 
thermal field, considering the complexity of their coupled dy- 
namics (cf. |[T3l ). and the anomalous response of the vortical 
condensate. 



Finally, we note that the cloud angular momentum (rota- 
tion rate) reaches its peak when the vortex leaves the conden- 
sate (f * 3200 eye), after which it undergoes a second near- 
exponential decay phase, the angular momentum of the field 
decaying such that by t w 5000 eye, only the thermodynamic 
fluctuations in (L z ) t h exhibited by the finite-temperature field 
at rest in the anisotropic potential remain. 



C. Heating of the atomic field 

We now consider the heating of the atomic field during the 
arrest of its rotation. As the system we evolve is Hamiltonian, 
with a time-independent potential, the total energy of the clas- 
sical field is a constant of the motion. Consequently, the trap 
anisotropy dissipates the angular momentum of the field by 
converting the rotational kinetic energy of the field into in- 
ternal energy lfT3l . and we therefore expect some heating of 
the field to occur, due to this redistribution of energy. We can 
estimate the heating of the field as follows: The rotational en- 
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ergy of the gas is initially CE TO t)i - ^i^i ~ 2-7 x I0 3 ha> r . We 
compare this to the initial thermal energy, which should be 
reasonably well estimated by the energy added to the ground 
vortex state in forming the initial thermal state (Sec. II Bl, 
Eth ~ 0.1£g ~ 7.6 x 10 3 hu> r . Assuming a linear relation- 
ship between thermal energy and temperature (which should 
be valid for the low temperatures we consider here, see ll35ll ), 
we might therefore expect an increase in the field temperature 
of ~ 30% during the arrest of the condensate rotation. 

In order to quantify the heating and its development during 
the field evolution, we estimate the (effective) thermodynamic 
parameters of the field (chemical potential fi and temperature 
T) using the procedure of [23 ,26]. As in [23], we allow for the 
differential rotation between the thermal cloud and the frame 
in which the classical-field cutoff is effected (the laboratory 
frame, in the present case). We perform the fit to the wing 
of the (time-averaged) radial distribution n{r) of the classical 
field, i.e., we fit over the radii range r e [r_ + 0.5ro, r tp ] where 
r_ marks the minimum of the Hartree-Fock effective potential 
V e ff - ( m /2)[a>^ - £2 t 2 h ]r 2 + 2U2Dti{r) ]36l , and r tp is the semi- 
classical turning-point of the low-energy space L |23l . We fit 
to radial densities of the classical field averaged over 10 eye 



periods, and assume the values of £2 t h calculated in Sec. IIIB 
The resulting estimates for li and T are presented in Fig.|4|a) 
We find that the temperature (dashed line) is approximately 
constant during the first ~ 1000 cycles of the field evolution, 
corresponding to the initial arrest of the thermal field compo- 
nent. The initial angular momentum of the thermal component 
is small, comprising some ~ 25% of the total angular momen- 
tum of the field, and so any heating of the field due to the re- 
distribution of the associated rotational energy is possibly too 
small to resolve above the uncertainty in the temperature esti- 
mates. Beginning at t « 1100 eye, corresponding to the start 
of the vortex-decay phase, the temperature exhibits a steady, 
approximately linear increase. We associate this increase with 
the conversion of rotational kinetic energy to thermal energy, 
due to the action of the trap anisotropy, and presumably also 
as a result of the scattering of excitations by the vortex which 
produces the frictional effect ifTOl . There is a final sharp rise 
in temperature at t « 3000 - 3200 trap cycles. This rise is 
perhaps due to the final decay of the vortex into excitations 
at the surface of the condensate IfTOll . We note however that 
the field, and in particular the condensate surface, undergoes 
strong fluctuations as the system passes through this transition 
to the vortex-free state, with, e.g., multiple (ghost) vortices 
present at the condensate surface at times. This can be viewed 
as the system essentially 'reversing' through the surface-wave 
instability arising from the relative motion of condensate and 
thermal cloud, in which vortices spontaneously grow from 
surface excitations ll26l l37l l38l . The thermodynamic param- 
eters may therefore be ill-defined during this period. After 
this period of strong surface fluctuations, the temperature lev- 
els off at T ss 2lhoL> r /k B , corresponding to heating of « 50% 
during the arrest. We note that the chemical potential [solid 
line in Fig. |4|a)] exhibits a slight downward trend beginning 
around t x 1 100 eye, falling by ~ 0Afia> r , though we note that 
this change is of the same order as the variation in estimates 
for ii obtained at late times t > 3500 eye. 
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FIG. 4: (Color online) Heating of the field. Evolution of (a) the 
effective chemical potential and temperature and (b) the condensate 
fraction. 



In Fig. Qb) we plot the condensate fraction f c = 
no/ Yjk>o n k obtained from the procedure outline in Sec. Ill A' 



The fluctuations in estimates of this quantity are large, as is 
expected given the short time scale (10 eye) over which the 
appropriate averages are taken, and the nonequilibrium nature 
of the field. However the condensate fraction estimates exhibit 
a clear downwards trend as time proceeds, and the condensate 
fraction drops by ~ 3% during the arrest. 



IV. DEPENDENCE ON TRAP ANISOTROPY 

We now consider the effect of varying the trap anisotropy 
e on the behavior of the classical-field trajectories. We intu- 
itively expect the rate at which the trap dissipates the angular 
momentum of the thermal cloud to depend strongly on the 
magnitude of the anisotropy, and a quantitative model for this 
dependence in the case of a classical (Boltzmann) gas was 
presented in II14I1 . By contrast, the efficiency with which the 
thermal cloud extracts angular momentum from the conden- 
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sate is dictated by the (longitudinal) mutual-friction coeffi- 
cient, which depends on the temperature and chemical po- 
tential of the field, in addition to its microscopic properties 
(4] [TO] [39] . We therefore expect, as discussed in the vortex- 
continuum analysis of |[T3l . to explore different regimes of 
relaxation dynamics as we vary the trap anisotropy and con- 
sequently the relative strengths of vortex-cloud and cloud-trap 
friction. In this section we consider the evolution of the clas- 
sical field in the presence of anisotropies of different magni- 
tudes, and observe the resulting differences in the behavior of 
the fields. 



A. Weak damping 

In Fig. |5ja) we plot the evolution of the angular momentum 
of the condensate and thermal cloud in the presence of a weak 
trap anisotropy (e = 0.005). We observe that after an ini- 
tial decline over the first w 10 4 trap cycles, the cloud angular 
momentum, despite exhibiting large fluctuations, maintains 
a reasonably steady mean value which increases slowly over 
the period t ~ [10 4 ,3 x 10 4 ] trap cycles. During this period 
the angular momentum of the condensate drops by approxi- 
mately a factor of 2, and large fluctuations in the distribution 
of angular momentum between the two components are visi- 
ble. This behavior is again apparent in the vortex precession 
frequency and cloud rotation rate presented in Fig. [5jb). This 
shows that the cloud exhibits a steady-state nonequilibrium 
behavior, i.e., its rotation rate remains somewhere between 
that of the precessing condensate mode, and that of the trap 
(£2 tI =0). In this regime its angular momentum remains ap- 
proximately constant, as the rate at which it gains angular mo- 
mentum from the decaying vortex matches the rate at which 
it loses angular momentum to the trap. We expect the angu- 
lar velocity at which this balance occurs, and indeed whether 
such a regime occurs at all, to depend strongly on the rela- 
tive strengths of the vortex-cloud and cloud-trap friction. We 
note that the angular velocity of the cloud appears to slowly 
increase over time, as the condensate loses angular momen- 
tum and its angular velocity increases, shifting the cloud ro- 
tation rate at which the angular momentum transfer rates are 
balanced. At t w 3.2 x 10 4 eye the vortex's precession accel- 
erates as it approaches the condensate boundary, and during 
the period t e [3.2 x 10 4 , 3.4 x 10 4 ] eye the adiabaticity of the 
vortex dissipation is lost as the field enters the critical regime 
associated with nucleation of the vortex at the surface 11231 . In 
this regime the condensate surface is unstable and as noted in 



Sec. IIIC multiple 'ghost' vortices may be present simulta- 
neously at the condensate boundary. A large amount of an- 
gular momentum is transferred non-adiabatically to the ther- 
mal cloud during this period of criticality, while small surface- 
mode oscillations of the condensate persist until t w 4 x 10 4 
eye. The angular momentum is subsequently dissipated from 
the thermal cloud over a period of ~ 3 x 10 4 eye, in an oscil- 
latory but approximately exponential decay phase. 



B. Strong damping 

We now turn our attention to a scenario of vortex arrest due 
to the presence of a strong trap anisotropy (e = 0.1). The dy- 
namics in this case are strongly nonequilibrium and the arrest 
of the field's rotation occurs on a shorter time scale than the 
cases already considered, and so here we form the density ma- 
trix Eq. |7]i by averaging classical-field samples over a shorter 
period of two trap cycles. In this case the angular momen- 
tum of the thermal cloud [solid line in Fig.[3Jc)] is rapidly de- 
pleted, and actually fluctuates below zero by t » 400 trap cy- 
cles. The condensate is slower to respond; during this period 
the condensate angular momentum decreases by ~ 7%. Sub- 
sequently, the cloud angular momentum fluctuates strongly, 
and rises to ~ 2000ft (close to its initial value), as the vortex 
is rapidly expelled from the condensate. The angular momen- 
tum is then dissipated from the cloud over a period < 200 
eye, and thereafter the angular momentum of the field fluctu- 
ates about zero. The rapid expulsion of the vortex, and strong 
fluctuations of the cloud rotation, are again visible in the cal- 
culated angular velocities [Fig. |5jd)] . 

The behavior of the cloud angular momentum in this tra- 
jectory suggests that the cloud is overcritically damped by the 
trap anisotropy [14|. Its response to the anisotropy appears 
unhindered by its coupling to the condensate via the vortex 
core, and it quickly yields the angular momentum it acquires 
from the vortex to the trap, despite its small angular velocity 
^th ^0.1. The relaxation process in this case is violently non- 
adiabatic, as evidenced by the strong fluctuations in the cloud 
angular momentum. Indeed during the period t ~ [600, 700] 
eye large, long-wavelength surface oscillations are visible in 
the field as the vortex precesses rapidly near the condensate 
boundary. It appears that the strong trap anisotropy and rapid 
'stirring' motion of the vortex conspire to strongly perturb the 
condensate bulk in this regime. During this period the mea- 
sured condensate fraction is suppressed [inset to Fig.|5jd)] due 
to the strongly nonequilibrium behavior of the condensate, in 
which surface-wave excitations define frames of rotation dis- 
tinct from that of the vortex. The decomposition of the field 
into condensed and noncondensed components must therefore 
be viewed with some caution in this strongly nonequilibrium 
scenario. Nevertheless, it is clear that the effect of the trap 
on the cloud dominates the vortex-cloud coupling in this sce- 
nario, in stark contrast to the near-adiabatic decay scenario of 
Sec. |IVA] 



C. Decay times 

We now consider the dependence of the relaxation times 
on the strength of the trap anisotropy. As the relaxation of 
both the condensate and the thermal cloud is generally nonex- 
ponential, we consider the times at which the condensate be- 
comes irrotational, and at which the total angular momentum 
of the field is lost. Due to the persistence of surface excitations 
which prevent the condensate angular momentum from reach- 
ing zero, we define the condensate stopping time t coi1c j as that 
at which the angular momentum of the condensate first drops 
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FIG. 5: (Color online) (a) Angular momenta and (b) angular velocities of field components in a simulation with trap anisotropy e = 0.005. 
(c) Angular momenta and (d) angular velocities of field components in a simulation with trap anisotropy e = 0.100. Inset: Evolution of the 
condensate fraction in the case e = 0.100. 



below 0.02H per particle, and the field stopping time as the first 
subsequent time Tfi e i<j at which the cloud angular momentum 
reaches zero. These times are presented in Fig.|6] We observe 
that the condensate arrest times (circles in Fig. |6t vary by two 
orders of magnitude over the range of trap anisotropics we 
consider. Moreover, the total-field relaxation times (plusses) 
become increasingly longer than the vortex relaxation times 
as the trap anisotropy is weakened, causing slower dissipa- 
tion of the cloud angular momentum even in the absence of 
the vortex. We therefore consider the time At = t^j - t coik j 
over which the angular momentum of the cloud dissipates fol- 
lowing the expulsion of the vortex. Although the angular mo- 
mentum lost in this final damping phase varies over the range 
(L.) t h ~ 3000-5000/z, precluding a precise analysis, the times 
At provide a useful characterization of the dependence of the 
cloud relaxation on the trap anisotropy. We perform a linear 
fit to At as a function of e in log-log space (inset to Fig. 
and find the scaling At oc e -2 - 1 , in good agreement with the 
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V. SUMMARY AND CONCLUSIONS 

We have carried out the first simulations of the arrest of 
a rotating Bose-Einstein condensate due to the presence of 
a trap anisotropy which includes the coupled nonequilibrium 
dynamics of the condensate and thermal cloud. Our method 
makes no assumptions of stationarity of a thermal bath, nor 
are the two components artificially given disparate rotational 
parameters. Rather, our approach describes the dynamical mi- 
gration of an equilibrium rotating thermal state to a new, irro- 
tational equilibrium, solely due to the action of the trapping 
potential. 

We observe for all parameters we considered that the rota- 
tion rates of the condensate and thermal component are dis- 
tinct during the decay. The anomalous rotational response 
of the precessing-vortex condensate can lead to a counter- 
intuitive spin-up of the thermal component during the decay, 
and we observe nonequilibrium oscillations in transfer of an- 
gular momentum between the two components as the vortex 
responds to the dissipative effect of the thermal cloud, which 
is itself damped by the trapping potential anisotropy. 

For trap anisotropies that are weak, the thermal field settles 
to a nonequilibrium steady state, with rotation rate interme- 
diate between that of the condensate and that of the (static) 
trap. In this scenario the angular momentum of the conden- 
sate is slowly depleted while that of the thermal cloud remains 
nearly constant, until the vortex nears the condensate bound- 
ary and the linearity of the vortex decay breaks down. For 



stronger trap anisotropies, the angular momentum of the ther- 
mal cloud may be almost entirely depleted before the conden- 
sate responds. 

We quantified the heating of the atomic field during the ar- 
rest, and found it to be commensurate with the conversion 
of rotational kinetic energy into thermal energy by the trap 
anisotropy. We also considered the time scales over which 
the condensate angular momentum and total field angular mo- 
mentum were dissipated, and found reasonable agreement be- 
tween the scaling of the thermal cloud spin-down time with 
trap anisotropy and the predictions of a Boltzmann gas model 

m. 

Possible extensions of this work include the description of 
vortex-lattice arrest, and the effects of temperature on the de- 
cay dynamics. In particular it would be interesting to study 
how the nonequilibrium rotational dynamics of the conden- 
sate approach those of a rigid body as the vortex density in- 
creases, and the effect of temperature on the cloud rotation in 
the steady-state scenario. 
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